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Abstract 

Motivated by the theoretical study of Saito and Ono (2011) on three crys- 
talhne forms of Si02 under pressure, quantum-chemical calculations on var- 
ious free-space clusters of SiO„ and GeO„ for 1 < n < 6 are reported here. 
Both neutral and cationic clusters have been examined, for both geometry 
and equilibrium bond lengths. Coupled clusters and correlation-corrected 
MP2 calculations are presented. For the cations, we emphasize especially 
the structural distortions occurring in removing degeneracies. 

FAGS: 31.15.Ne, 36.40. Qv 
1. Introduction 

The background to the present study is to be found in the theoretical 
work of Stenhouse et al. in this Journal |JL]. This dealt with scattering 
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Figure 1: (Color online) Shows Si atom approaching O2 molecule. 



intensities and partial structure factors in vitreous silica, both X-ray and 
neutron diffraction experiments being invoked. Our interest in the present 
Letter focusses on the electron density distribution in free-space clusters such 
as Si02, Si04, and SiOe- This quantum-chemical study has been directly 
motivated by the theoretical work of Saito and Ono 2'1 on different crystalline 



forms of Si02 under pressure (see also [Sj on vitreous silica). In addition to 
our calculations on the above three neutral clusters of the form SiO„, we 
have also considered cationic free-space clusters of this kind. 

2. Free-space cluster Si02 by Hartree-Fock (HF) theory 

All the clusters studied in this work have been first considered in free 
space, by performing geometry optimization within the coupled cluster ap- 
proximation with single and double excitations (CCSD) 0,|5[|6[. Correlation- 
corrected MP2 calculations were then carried out for the optimization. In all 
the calculations, performed by means of the Gaussian 09 software [3], the 6- 
31 basis set, supplemented by polarization and diffuse functions (6-31G+**), 
was used. 

Fig. [1] shows the way in which an O2 molecule, at a fixed bond length 
1.21 A, is brought up to a single Si atom in the neutral cluster Si02. The 
distance d was varied from 2 A to 1.75 A at intervals of 0.05 A. Table [1] 
reports the HF energies for singlet, triplet, and quintet spin states. While 
O2 still retains its spin density at the largest distance di = 2 A, the triplet 
state energy lying below the singlet value, a crossover in energy taking place 
around 1.85 A; at 1.75 A separation the singlet state clearly lies lowest. 
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di (A) 




E (Hartree) 






singlet 


triplet 


quintet 


2.00 


-439.104 


-439.118 


-439.013 


1.95 


-439.112 


-439.122 


-439.002 


1.90 


-439.121 


-439.125 


-438.990 


1.85 


-439.129 


-439.127 


-438.976 


1.80 


-439.137 


-439.128 


-438.960 


1.75 


-439.144 


-439.128 


-438.941 



Table 1: Energies in Fig. [T] as a function of distance c?i, for different spin states. 







distance (A) 








angle (°) 


1- 


-2 


1.729 


2- 


-1- 


-3 


82.25 


1- 


3 


1.729 


2- 


-1- 


-4 


108.5 


1- 


-4 


1.639 


2- 


-1- 


-5 


108.5 


1- 


-5 


1.639 


4- 


-1- 


-5 


130.14 



Table 2: Bond lengths and angles for optimized low- lying isomer of Si04, as shown in 
Fig. El 

3. Neutral free-space clusters Si04 

Fig. [21 calculated by HF theory, shows what we consider to be the low- 
lying isomer of the free-space cluster Si04. The geometry of this present 
proposal for the ground state is recorded in Table [2j It is seen that there 
are two SiO distances which we determine to be 1.64 A and 1.73 A respec- 
tively. The corresponding angles range from 82.3° to 130°, via the (almost) 
tetrahedral angle of 109°. It is already relevant to make some contact with 
the crystalline calculations of Saito and Ono ^ (see especially their Ta- 
ble II). These workers show that at ambient temperature and pressure, the 
ground-state structure a-quartz (g-Si02 in their Table II) has two slightly 
different bond lengths 1.60 A and 1.61 A, with angles near the tetrahedral 
value ~ 110°. Above a pressure of 2 GPa, a rutile structure (r-Si02) is 
formed, and Table II records two SiO bond lengths of 1.75 A and 1.79 A, 
which are at least fairly near to our predicted second bond length 1.73 A for 
the free-space Si04 cluster in Fig. O Angles were not recorded however for 
r-Si02 in their Table II. 

We shall return briefly to our further results pertaining to Si04 clusters 
when we summarize our calculations on cationic free-space clusters below. 
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Figure 2: (Color online) Si04 optimized low-lying isomer. For precise bond lengths and 
angles, see Tabled 




Figure 3: (Color online) Shows two different views of isomer of Si04, with geometry as in 
Table H 

4. Neutral free space cluster of a Si atom with six surrounding O 
atoms 

Fig. shows first a neutral cluster in which atom 1 is Si and tetrahedral 
angles for 4 surrounding O atoms are assumed. The energy and bond dis- 
tances are recorded in Table [3l the larger distance 1-2 and 1-3 being around 
2.3 ± 0.1 A, compared with 4 shorter bond lengths around 1.66 ± 0.01 A. 
Table H] shows how variation of the geometry away from the tetrahedral an- 
gle 109.5° entering Table [3] lowers the HF energy by ~ 0.008 Hartree, with 
the cluster form depicted now in Fig. Hb- We note briefiy that in Table II 
of Ref. [2I, the longest bond length recorded is for r-Si02, and is 1.81 A. We 
have also considered the stability of the neutral cluster SiOg with respect to 
an isolated Si atom plus two ozone molecules O3, finding a relative energy 
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(a) (b) 
Figure 4: (Color online) Proposed low-lying isomers of SiOg. 



n 


E{0) 


E{Si) 


^(Si + nO) 


E{SiOn) 


AE 


2 


-74.940 


-289.017 


-438.897 


-439.244 


-0.346 


4 


-74.940 


-289.017 


-588.777 


-589.384 


-0.607 



Table 3: Stability with respect to the elements of isomers Si02 and Si04 displayed in 
Figs. [T] and m respectively. Energies are in Hartrees. 

difference of AE = —0.373 hartrees. 

5. Some less detailed results with the Si atom replaced by Ge in 
free-space neutral clusters 

Prompted by the solid-state theoretical calculations presented by Saito 
and Ono j2|, we have also studied GeO„ clusters, though in somewhat less 
detail than their Si counterparts. One of the configurations of the free-space 
neutral cluster is depicted in Fig. |5]for GeOg. The bond distances and angles 
are collected in Table [5l for comparison with SiOe for the neutral cluster 
shown in Fig. Hb- 

6. Singly-charged positive ion clusters SiO^ and GeO^ 

Having treated the neutral clusters of SiO„ in some detail, we felt it of 
interest in longer terms to estimate bond length distortions due to ionization 
for both Si04 and Ge04 . We note first that if we were dealing with the 
free-space tetrahedral molecule silane, then we could appeal to experiment 
for the symmetry of the cation SiH4 (see, e.g., Ref. j^). 

We show therefore in Fig.[6]the proposed structures for these neutral Si04 
and cationic clusters Si04 from the present quantum-chemical technique. 
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distance (A) 




n Ti crl fi ( 1 


SlOe 


(Fi} 


?. m, 


E = 




739.4260 Hartree 


1-4 




1.647 




5- 


-1-4 109.5 


1-5 




1.669 




5- 


-1-6 109.5 


1-6 




1.67 




5 


-1-7 109.5 


1-7 




1.679 








1-2 




2.195 








1-3 




2.377 








SiOe 




?.llb), 


E = 




739.4343 Hartree 


1-4 




1.664 




5- 


-1-4 123.07 


1-5 




1.664 




5- 


-1-6 84.77 


1-6 




1.664 




5- 


-1-7 123.07 


1-7 




1.664 








1-2 




2.208 








1-3 




2.208 









Table 4: Shows geometries of low-lying isomers of SiOe depicted in Fig. 21 




Figure 5: (Color online) Low-lying isomer for GeOg, with bond lengths and angles as in 
Tab. El 

Bond lengths and angles are correspondingly recorded in Tab. [61 Analogous 
results are then shown in Fig. [7] and Tab. [7] for neutral Ge04 and cationic 
clusters GeO^. In the light of the earlier discussion, we note in particular 
that now two bond lengths are in evidence, the longer being 2.2 A for GeO^. 

Additionally in Table[8]we record our calculations of the Mulliken charges, 
which demonstrate the expected shift of electronic charge from Ge to the O 
nuclei, due to the high electronegativity of the latter atom. 

We conclude this section by recording for both Ge04 and Ge04 the one- 
electron self-consistent field eigenvalues for these neutral and also charged 
clusters (Tab. [9]). As with the bond lengths, even though we are not dealing 
with stable gas phase molecules like SiH4 and its cation, there are parallels 



6 



distance (A) 



angle (°) 



1- 7 

2- 7 

3- 7 

4- 7 

5- 7 

6- 7 



2.333 
2.333 
1.797 
1.797 
1.796 
1.797 



3-7-4 
3-7-5 
3-7-6 



123.54 
130.04 
78.84 



Table 5: Bond lengths and angles for optimized low-lying isomer of GeOg, as shown in 



Figure 6: (Color online) Low-lying isomers predicted theoretically for neutral Si04 
and cationic Si04 , the total ground-state energies being E — —589.3844 and E = 
—588.8507 a.u., respectively. Bond lengths and angles are given in Tab. HI 

with the Jahn- Teller removal of degeneracy. 

7. Summary 

Our major conclusion from the present theoretical study concerns the 
predicted low- lying isomers for the neutral free-space clusters SiO„ (1 < ??, < 
6), and in briefer fashion for the corresponding Ge clusters. As stressed 
in Sec. [H the immediate motivation was afforded by the recent theoretical 
work of Saito and Ono [2] on three different crystal structures of Si02 under 
pressure (see also our own early study on diffraction from vitreous silica in 
p|). It is to be hoped that our present results will stimulate experimental 
studies on the free-space clusters SiO„ for small n, and possibly also on 
the comparison between the electron density distribution in crystalline Si02 
under pressure with vitreous Si02, also under compression. 



Fig. El 




7 



distance (A) angle (°) 



Si04 

1-2 
13 

1-4 
1-5 



1.63 
1.63 
1.63 
1.63 



2-1-4 
2-1-5 
2-1-3 



137.12 
137.13 
62.24 



SiO+ 

1-2 

13 

1-4 

1-5 



1.729 
1.729 
1.639 
1.639 



2-1-4 
2-1-5 
2-1-3 



108.5 
108.52 
82.25 



Table 6: Bond lengths and angles for neutral cluster Si04 and cation Si04 , as shown in 
Fig. El 
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distance (A) 








angle (°) 


Ge04 












1-2 


1.770 


2- 


-1- 


-4 


140.31 


13 


1.770 


2- 


-1- 


-5 


57.32 


1-4 


1.770 


2- 


-1- 


-3 


140.35 


1-5 


1.770 










GeO+ 












1-2 


1.755 


2- 


-1- 


-4 


120.58 


1-3 


1.755 


2- 


-1- 


-5 


95.9 


1-4 


1.755 


2- 


-1- 


-3 


118.13 


1-5 


2.252 











Table 7: Bond lengths and angles for neutral cluster Ge04 and cation Ge04 , as shown in 

Fig.m 





Si04 




SiO+ 




SiOe 




Ge04 




GeO+ 


Si 


2.458 


Si 


2.202 


Si 


2.440 


Ge 


2.416 


Ge 


2.348 





-0.615 





-0.126 





-0.600 





-0.604 





-0.420 





-0.615 





-0.126 





-0.600 





-0.604 





-0.421 





-0.615 





-0.474 





-0.020 





-0.604 





-0.424 





-0.615 





-0.474 





-0.600 





-0.604 





-0.084 



O -0.600 
O -0.020 



Table 8: Mulliken charges in free space neutral and cationic clusters. 
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GcO 


4 


GeOt 


-405.39349 


-1.79674 


-405.63174 


-20.88425 


-2.03094 


-0 


95041 


-52.29426 


-1.79585 


-405.63137 


-20.88342 


-2.03063 


-0 


93999 


-46.38488 


-1.44388 


-52.53368 


-7.58331 


-2.02971 


-0 


93207 


-46.38488 


-1.43636 


-52.53365 


-7.58274 


-2.02948 


-0 


90334 


-46.38119 


-1.2333 


-46.62474 


-5.56251 


-1.6546 


-0 


89927 


-20.68764 


-1.23317 


-46.62464 


-5.56214 


-1.64433 


-0 


86177 


-20.68758 


-0.7747 


-46.62181 


-5.55912 


-1.63354 


-0 


84468 


-20.6875 


-0.67349 


-46.62164 


-5.55908 


-1.57882 


-0 


83861 


-20.68743 


-0.65732 


-46.62148 


-5.55805 


-1.57679 


-0 


83284 


-7.34623 


-0.65728 


-46.62147 


-5.5578 


-1.48633 


-0 


83033 


-5.32489 


-0.56935 


-21.00925 


-2.03818 


-1.47287 


-0 


81215 


-5.32489 


-0.56444 


-21.00894 


-2.03762 


-1.4617 


-0 


80145 


-5.32224 


-0.55312 


-20.92582 


-2.03492 


-0.9934 


-0 


78798 


-1.80179 


-0.55302 


-20.92566 


-2.03362 


-0.98829 


-0 


78047 


-1.80075 


-0.53758 


-20.9248 


-2.03347 


-0.98721 


-0 


77532 


-1.80075 


-0.52855 


-20.88442 


-2.03286 


-0.9819 







Table 9: One-electron eigenvalues of both Ge04 and Ge04 (atomic units). 



11 



